Interstellar Extinction in the Vicinity of the Galactic Center 
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ABSTRACT 

We present J (1.2 pm), H (1.6 pm), K' (2.11 pm) and Br 7 (2.166 pm) 
images from four large regions within the central 40 pc of the Galaxy. Localized 
variations in the extinction, as determined by observations of the stellar 
population, are examined using the median H-K' color as a function of position 
within each region. As the value of the derived extinction from the stars is 
dependent upon the intrinsic magnitude of the assumed stellar type, the J-H 
vs. H-K' diagrams are first used to investigate the distribution of stellar types 
in the four regions. We find that there is a distinct OB population, contrary 
to earlier assumptions, with the ratio of K and M giants and supergiants to 
OB stars approximately twice that of the solar neighborhood. Although the 
on the scale of ~1' fluctuations in the extinction are on the order of Ay ~2, 
throughout the entire region the extinction varies from Ay^, 25 to Ay ^ 40. 
We also examine whether there is any variation in the extinction and stellar 
number density relative to the unsual radio features in these regions and do not 
find a significant correlation. 

Spectral imaging in Br 7 2.166 pm emission shows a strong morphological 
correspondence between the 6 cm radio images and the diffuse Br 7 emission. 
By comparing the theoretical Br 7 flux derived from the radio flux using 
recombination theory, with our measured Br 7 flux, we obtain a second, 
independent, estimate of the extinction. We compare the two data sets and find 
that the extinction as derived from the stars is consistently greater, sometimes 
by a factor of two, than the value of the extinction derived from the Bry images. 
The derived extinction in various regions, however, is insufficient for any of these 
regions to be located behind the molecular clouds as previously observed in the 
Galactic Center, consistent with the theory that the observed radio emission is 
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produced on the foreground surface of these clouds. 


Subject headings: Galaxy: center-interstellar extinction 
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1. Introduction 

The Galactic Center (GC) is unobservable at visible wavelengths because of interstellar 
extinction, usually assumed to be A v ~ 30, but known to vary between 5^^ ^ 60 
(Catchpole, Whitelock & Glass 1990). The bulk of the extinction is believed to be the 
result of clouds of dust and gas within 4 kpc of earth and not due to clouds within the 
inner few hundred parsecs of the Galaxy (Sanders, Scoville & Solomon 1985). Previous 
work to quantify the extinction towards the GC went to a limiting magnitude of K=12.1, 
H=12.8 and J=13.5 (Catchpole et al. 1990), which, as shall be shown in §5, introduces a 
bias in the derived extinction. Other photometric studies in the GC have gone to fainter 
magnitudes, but to date these studies have been limited to the inner few parsecs (e.g. Blum 
et al. 1996). This work therefore removes the sample bias of previous studies and provides 
an examination of the variation in extinction and stellar population within ^5' of Sgr A. 

In addition, this study examines whether or not there is a correlation between variations 
in the extinction and the position of many of the unusual radio continuum features which 
make the inner ~40 pc of the Galaxy so intriguing. 

Prior determinations of the large scale extinction variations derived from stellar colors 
have begun with the assumption that the stars observed in the GC are predominately 
evolved K and M giants. The validity of this assumption has been questioned for some time, 
at least since Lebofsky, Rieke &; Tokunaga (1982) first suggested that they were observing 
M supergiants in the GC and proposed there had been a recent burst of star formation. 
Several of their M supergiant candidates were subsequently shown to be M giants (Sellgren 
et al. 1987), but IRS 7 was confirmed to be a young M supergiant (Rieke & Rieke 1989). 
Therefore the assumption deemed most appropriate was one which included a population 
of M giants and supergiants. That this is not the case, however, was suspected when young 
Of star candidates were identified within the central parsec (e.g. Forrest et al. 1989; Allen, 
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Hyland & Hillier 1990; Krabbe et al. 1991). When additional young massive WN and 
OB stars outside of the central parsec were unequivocally identified by their spectra (e.g. 
Cotera et al. 1996 (Paper I); Figer, McLean & Morris 1996), the assumption of a population 
of exclusively K and M giants and supergiants was no longer valid. Recently, larger scale 
NIR observations of the central parsecs have confirmed that there is a luminous young 
population throughout the central 30 pcs (Mezger et al. 1999, Philipp et al. 1999) 

With the discoveries of O, B and WN star candidates at several locations within the 
central ~40 pc (e.g. Cotera et al. 1994; Nagata et al. 1995; Figer, McLean & Morris 1995; 
Blum, DePoy & Sellgren 1995; Paper I; Cotera et al. 1999 (Paper II)), in addition to a 
renewed interest in the general character of the stellar populations in the GC (Mezger et al. 
1999, Philipp et al. 1999), there has also been a renewed effort to determine the intrinsic 
stellar properties of these stars (e.g. Najarro et al. 1997). In order to accurately address 
these and other areas of GC research, a good understanding of the extinction is necessary. 
We have a degenerate problem however: in order to use the infrared color-magnitude 
diagrams to determine the extinction, we must assume an intrinsic stellar type- but which 
stellar type to use is no longer obvious. Therefore, before we can determine the extinction 
from the H-K' colors, in §3.2 we use color-color diagrams to estimate a realistic ratio of K 
and M giants and supergiants to OB stars, which we then employ in our estimate of the 
extinction. 

In addition to the photometric data, we have used the difference in the observed Br 7 
(2.166 £tm) emission from several radio continuum regions, to that predicted from Case 
B recombination theory to derive a value of the extinction specifically for these features. 
Comparison shows that the extinction derived from radio recombination is up to a factor of 
2 less than that derived from the stellar photometry. The Br 7 emission images with their 
resultant extinction estimates are presented in §4. In §5 we discuss the implications of the 
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two data sets taken together and argue that the radio emission must be from the foreground 
surface of the molecular clouds in the region and the stars must also be in front of these 
clouds. 


2. Observations and Data Reduction 

Observations were made at the Anglo-Australian Observatory on 1992 July 13-14 and 
1993 June 13-15 with the facility Infrared Imaging Spectrometer (IRIS) which utilizes a 
Rockwell 128x128 pixel HgCdTe (NICMOS II) array. The 0.6"/pixel scale was selected 
for all observations. IRIS uses the K' filter (central wavelength 2.1 /rm; Wainscoat & 
Cowie 1992), and standard H (1.6 gm) and J (1.2 /rm) filters for broad band observations 
(differences in the photometric measurements of K and K' are briefly discussed in §3.1). All 
data were linearized, bias subtracted, flat fielded and sky-subtracted with the astronomical 
software reduction package FIGARO. 

We obtained J, H and K' images in four regions shown as shaded areas in Figure 1. 
The four regions are: 1) a portion of the Straight Filaments including AFGL 2004 (the 
‘Quintuplet cluster’), GO. 15-0. 04 (the ‘Pistol ’) and GO. 18-0.04 (the ‘Sickle’, Yusef-Zadeh & 
Morris 1987b), hereafter referred to as QPS; 2 ) the Arched Filaments including the El & E2 
Filaments and GO. 10+0.02, hereafter referred to as the El & E2 Filaments, (nomenclature 
from Morris & Yusef-Zadeh 1989); 3 ) the region south of the Arched Filaments surrounding 
the radio emission regions H1-H8 (nomenclature Yusef-Zadeh &: Morris 1987a); and 4 ) the 
radio emission regions A-D just east of Sgr A East (Yusef-Zadeh & Morris 1987a), which 
Zhao et al. (1993) call the Sgr A East H II cluster, hereafter referred to as Sgr A East. 

A combined broad band J, H and K' image of the El &; E2 Filaments was presented as 
Figure 1 of Paper I, the remaining three images are presented as Figures 2-4 overlaid with 
the available 6 cm radio contour maps. The figures are mosaics of 70" x 70" J, H, and K' 
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images which are combined in a three-color composite. These images were obtained during 
generally good conditions, with the image of QPS obtained on a night of £0.8" seeing. For 
the other images, the seeing was ~1.0-1.5". Astrometry was performed using SAO stars 
within the large images. In most of the images, however, only two SAO stars were available 
which did not allow us to simultaneously constrain the three free parameters: RA plate 
scale, Dec plate scale and rotation angle. Therefore, on the QPS, El and E2 Filaments, 
and H1-H8 images, positions are accurate to £2.5". On the image of Sgr A East, optical 
astrometry was available for enough of the stars that the positions are good to ~ 1". 

Photometry results for all of the images were obtained using the crowded field 
procedures suggested for use with IRAF/DAOPHOT. Photometry was performed on the 
individual frames rather than the large images, to improve accuracy and provide a means 
of determining the errors in the derived magnitudes from the overlapping regions. Each 
individual 128x128 pixel image is extremely crowded; for instance, in each of the 70" x 70" 
images at K' there are £400 stars. This crowding, the result of the large number of stars 
we observe while looking through the galactic plane, is the primary cause of the derived 
errors in magnitudes. Photometric errors are discussed in detail in §3.1. Only one frame 
per wavelength was taken at each position, so only one measurement is available for most 
stars, which contributes to the uncertainty: the magnitude error is thus 0.16 in K' for the 
worst case, although the typical errors are £0.1. 

Our sensitivity for these observations varies slightly for each region as the regions were 
taken on different nights and in some cases different years. Plots of the log of the number of 
stars as a function of magnitude for each of the four regions are presented in Figure 5. Since 
we are dealing with a combination of disk and bulge stars, we know from Cohen (1993), that 
in the magnitude range we are observing, the log N vs. magnitude is a relatively smoothly 
increasing function with decreasing magnitude. Therefore we fit a second order polynomial 
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to the regions of the plot unaffected by either saturation or detection limits: in essence 
the well sampled regions of the plot. In order to determine the limiting magnitudes first 
we determined the standard deviation (a) of the fit in the well sampled region . We then 
used the point at which the data deviated from the fit by ^ 5a as the limiting magnitude 
in H and K' and ^ 3a for J as the fewer number of stars detected in J increased the noise 
in the fit. The values for J, H and K' determined in this way for each of the regions are 
presented in Table 1. The determination of the limiting magnitude at K' in Sgr A East is 
effected by the significant noise in the region used for the fit due to the smaller sample, 
therefore the limiting K' magnitude presented in Table 1 is taken to be 3a. In addition to 
the limiting magnitudes, Figure 5 illustrates the variations in the number of stars detected 
in the different regions. The El & E2 Filaments and H1-H8 have approximately the same 
numbers of stars, which follows as they cover similar areas on the sky, and were taken in 
similar conditions. The data set around Sgr A East has considerably fewer stars due to the 
fact that only a ~3'x3' region is included as opposed to ~4'x8' in the El & E2 Filaments. 
There are significantly more stars detected in QPS due to the excellent seeing during the 
night of 14 July 1992 and notably lower limiting magnitudes in H and K'. 

Spectroscopic observations of the Pistol, a part of the Sickle, GO. 10+0.02, part of 
the E2 Filament, HI, H2, H5-H7, H8, and Sgr A East regions A-D utilized data cubes, 
with axes a, 5 and A. The data cubes were obtained by observing in spectral mode with 
a 70 ,/ x2.0" slit (resolution R~250) at fixed right ascension while the telescope stepped in 
declination at a rate of 36 "/hr, producing pixels which are approximately square in a 
and 5 with a 0.6"/pixel plate scale. The cubes were flat fielded, straightened and had the 
atmosphere removed by dividing by a main sequence G standard, generally taken before 
and after each data cube was obtained. Spectra of main-sequence early G standards were 
used for all flux calibrations (Allen & Craig 1983). The flux calibration routine assumes the 
star radiates as a black body; a temperature of 5600 K was used for all G standards (Allen, 
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D. A. 1976). Spectral images of each source were then extracted from the data cubes by 
coadding the desired a, 5 planes (usually 3 planes centered on the emission wavelength) and 
using adjacent planes, containing no emission line features, for sky continuum subtraction. 

Since the steps in declination used to build up the data cubes were not precisely 0.60", 
a result of the limits of telescope pointing accuracy, and there was some drift in RA over the 
often hour long integrations, the images extracted from the cubes are slightly irregular in 
RA and Dec. Therefore, the radio map overlays used in all figures do not always correspond 
precisely to the Bry emission images. In addition, since the astrometry of the stellar images 
is only good to ~ 2.5", and there is a very strong morphological correspondence between 
the radio continuum and the Br 7 emission, the images which include Br 7 emission use the 
correspondence between the radio continuum and the Br 7 emission to establish positions 
which then are accurate to <;1". Continuum images were also extracted from the cubes to 
determine the position of stars relative to the Br 7 emission and the broad band H and K' 
images were then overlaid using these stellar positions. Due to the irregular pixel size in the 
cubes, the positions of the stars as seen in the continuum images extracted from the cube 
do not always correspond exactly with the broad band H and K' images. 


3. Photometry 
3.1. K vs. K' Errors 

IRIS utilizes the K' filter (2.11 fim) rather than the traditional K filter (2.2 /mi). The 
K' band as described by Wainscoat & Cowie (1992), was selected because it reduces the 
thermal component of the background allowing for deeper imaging in the same integration 
time. In Wainscoat & Cowie (1992), they derive both empirical and observational relations 
for the conversion from K' to K, of the form 
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K' - K — a(H - K ) (1) 

where a needs to be determined for out specific observations. Wainscoat & Cowie (1992) 
derive a value for a of 0.19 using the central wavelengths of the different bands, and 
0.22±0.04 based on observations of nearby stars; when they include two highly reddened 
stars a — 0.18 ± 0.04. For our case, although additional observational calibration would be 
best, the value which includes the reddened stars is the appropriate choice. One additional 
problem, however, is the fact that the observationally calibrated value was based on 
observations performed at Mauna Kea while our observations were taken at the AAT. 

In order to assess the effect on the derived relation due to the difference in observatory 
locations, we have used the atmospheric transmission code ATRAN (Lord 1992) to simulate 
the atmospheric conditions at the two sites. We have convolved these atmospheric models 
with the filter functions for H, K and K 7 . Additionally, in order to investigate any differences 
which might arise from stellar type and extinction, we have further convolved the models 
with hot and cool blackbody spectra and the extinction curve of He et al. (1995). We 
find no theoretical dependence based on stellar type or extinction, but there is a correct 
required due to atmospheric differences: going from the AAT and Mauna Kea we find 
<5a: = 0.013 ± 0.001. Therefore, adjusting the most appropriate value from Wainscoat & 
Cowie (1992) to the AAT, we use c*=0.17±0.04 as the most appropriate value, throughout 
this paper. We note, however, that absolute values for a derived in the analysis based on 
the blackbody curves, were lower than 0.17, which indicates that the values of K', and 
quantities derived using this conversion may be overestimates. The value used throughout 
this paper for a represents what we believe to be the best estimate possible, given the 
various constraints. 


In determining the extinction to our sources, we use the extinction law of He et al. 
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(1995) which we analytically transpose to K' using equation (1) and our derived value for 
a. We adopt the value E^-v/E B -v = -2.72 throughout this paper. 

Although the DAOPHOT routine assigns errors to each of the stars based on the 
variations in the background and fit to the PSF determined for each individual frame from 
stars in that frame, in regions of high confusion such as the GC, it is advisable to make 
an independent determination of the errors by comparing the values derived for the same 
stars in different frames. We note that the largest errors are seen in the QPS data. On 
the night that the QPS data was taken, the seeing was excellent, <;0.8". Since our pixel 
size is ~0.6"/pixel, the data for QPS were undersampled, which results in larger errors. 
The difficulty of DAOPHOT in fitting stars with a PSF taken from undersampled data is 
reflected in errors derived by DAOPHOT which are approximately twice as large for the 
QPS stars as for the stars in all the other regions. Each of the images (Figures 2-4) consists 
of a mosaic of 70"x70" images. These individual frames overlap the adjacent images along 
their edges by ~18" (~30 pixels). We therefore used stars detected in two images, as well 
as those in the corners where four frames overlapped, to obtain additional estimates of the 
photometry errors. One caveat to this approach is that as we transition from examining all 
stars detected, to looking only at stars found in two overlapping frames, to using only stars 
found in four overlapping frames, fewer and fewer stars are examined such that at some 
point low numbers statistics must be considered. In the large mosaic images (QPS, El and 
E2 Filaments, and H1-H8) there were ^1000 stars detected in both H and K' in regions 
where two frames overlapped and ^300 stars found in J, H, and K 7 for the same regions. In 
regions where four frames overlapped, however, ~65-100 stars were detected in H and K 7 , 
and only 13-23 for stars found in J, H and K'. 

The median error as a function of magnitude for all of the images, using the errors 
determined from stars found in two overlapping frames, is shown in Figure 6. As expected, 
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the error increases as we go to fainter magnitudes. The median error was selected in the 
following analysis, as opposed to the mean, to reduce the influence of random outliers. 

In addition to taking 2 or 4 overlapping frames to determine the errors, three separate 
cases, labeled in column 1 of Tables 2a and 2b, were investigated according to the data 
to be presented. A: stars identified in both the H and K 7 bands. B: stars found in both 
the J and K' bands. C: stars which appear at all three wavelengths. The median error 
determined for the four large images for various cases are presented in Tables 2a and 2b. 
As the table shows, limiting the stars used for analysis to those found in J, H, and K 7 , 
generally improves the errors for the stars in H and K' as the selection criteria eliminates 
the fainter stars which typically have larger errors. This is consist with the results presented 
in Figure 6. Also notable in Table 2a, the error determined using four measurements of the 
same star (cases 4A-4C) is ~0. 01-0. 06 magnitudes larger than those determined from only 
two overlapping frames. This could be a result of the much smaller sample size available 
to determine the median error in the case of stars with four measured magnitudes (e.g., in 
the QPS there are 1659 stars in case 2A and only 78 in case 4A). The median error for the 
same stars as determined by DAOPHOT are also given in Tables 2a and 2b. The errors 
determined by DAOPHOT are consistently lower than those determined by comparing the 
same stars in different frames. In all subsequent analysis we use the errors as determined 
from two overlapping frames for the appropriate situation (e.g. case C for the J-H vs. H-K' 
plot of QPS), as they provide a more conservative error than those derived by DAOPHOT, 
and do not suffer from the small number statistics such as those determined from four 
overlapping frames. The specific case used to establish the error bars is noted on all plots 
and corresponds to the values in Tables 2a and 2b. 
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3.2. Color-Color Diagrams 

The fundamental problem of using broad band colors to determine the interstellar 
extinction in regions of high extinction is having two independent color parameters to 
constrain: the intrinsic colors of the stars, and the color due to interstellar extinction. 
Previously (e.g. Catchpole et al. 1990), the assumption has been that since the population 
of stars in the Bulge are old stars, no more OB stars can exist in the GC and thus the 
population is assumed to be mixture of K and M giant and supergiants (which we will call 
simply KM stars). Recently however, NIR spectroscopy of GC stars has demonstrated that 
there are hot young stars within the central ~50 pc of the Galaxy (e.g. Krabbe et al. 1991, 
Cotera et al. 1994; Figer et al. 1995; Blum, et al. 1995, Paper I, Paper II). Therefore we 
must use the available data to examine the assumption of a strictly KM star population 
and attempt to derive a more realistic determination of actual stellar population. 

We begin to constrain the intrinsic color of stellar population at the Galactic Center by 
noting that the only stellar types which can be detected at a distance of 8.0 kpc (Reid 1993), 
assuming a minumum of Ay=20 and using our observational limits, are stars with intrinsic 
M j <;-2.0. For reference, the intrinsic magnitudes of numerous stellar types are presented in 
Table 3. Using the derived limit, we immediately see that we are observationally restricted 
to most types of OB stars (including those on the main sequence), all M giants and all K 
and M supergiants. However, the spectral types of K giants which can be detected in the 
J bands drops sharply as a function of increasing Ay- For Ay^l8, we can observe all K 
Giants to K0, at Ay~20, we can no longer detect spectral type K0, and at Ay~25, we can 
no longer detect any K giants. Therefore, if we use the J band data in determining our 
extinction values, our results would be skewed towards foreground stars. Therefore, we 
will use only the H-K' color to determine the interstellar extinction. As we noted above, 
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however, we have two color parameters to determine and it is impossible to constrain two 
parameters with only one available color. We are therefore forced to use the J band data, 
in conjunction with the H and K' data, to explore the relative stellar populations. By using 
the J band data, we are limiting the completeness of our sample; however , the assumption 
implied by this approach - the stellar population as measured by observations in J, H and 
K' is approximately equivalent to the actual population - is much better than the previous 
assumption of no OB stars. Therefore, although there are limitations to the method, it is 
the best possible solution given the constraints. 

The J-H vs. H-K' color-color diagrams for the four regions are presented in Figures 
7-8. On each plot we have included the extinction curve for KM and OB stars, using the 
extinction curve of He et al. (1995) where E//_y/E B _ v — —2.76, Ej_y /Eg-v — —2.25 and 
Ek'-v/Eb-v = —2.72 as explained above. Using the appropriate errors for the specific 
data set, we have also plotted 1 a error bars for the extinction curves, which are shown 
as shaded regions (see Tables 2a-2b for exact error values). The intrinsic magnitudes 
used as the initial point for the two lines are determined by using the average of intrinsic 
colors for the two different population groups based on Table 3, with Mjf converted to 
Mtf/ using eguation (1). Throughout the remainder of this paper for the intrinsic colors 
of OB and KM stars we use: (J r — Hi) OB = —0.097 ± 0.033, (J* — H iKM = 0.71 ± 0.15, 
(Hi - Kfi OB = -0.044 ± 0.025, {H t - K' i ) KM = 0.161 ± 0.037. 

As can be seen in Figures 7 and 8, many of the stars lie along the KM curve; however, 
in all the images except Sgr A East, a substantial number of stars also lie along the OB 
curve. In order to quantify these qualitative results, we use the extinction curves of KM 
and OB stars to determine which star are contained within the various shaded regions in 
Figures 7-8 (the extinction curve plus the estimated la errors). Based on the location of 
a star within these error boxes, a spectral type of either KM or OB is assigned. As can 
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be seen in Figure 7, in the QPS region some stars are located within the error box of both 
spectral types. No effort has been made to force these stars into one type or another, and 
the extent to which the samples overlap suggests that the results of the analysis should be 
treated with caution. Due to the incompleteness of the sample due to limits at J, these 
results should be treated as indicating trends only. Stars outside of the error boxes are 
assumed to have either infrared excess or blue excess and are eliminated from the sample. 

In order to assess whether or not the scatter in the J-H vs. H-K' plots is actually 
the result of sampling two distinct stellar populations or is simply the result of random 
photometric errors greater than the estimated values, we looked at the number of stars with 
infrared excess versus blue excess in relation to a strictly KM population and a two group 
KM and OB population. If the scatter we are seeing is due only to random errors, we would 
expect that there should be equal numbers falling outside the error boxes of the different 
populations. Looking at stars outside of the error box for the KM extinction curve only, 
we find there are 3. 5-5.0 times as many stars with red excess as those with blue excess. If, 
however, we divide the population into KM stars and OB stars, and look at the number of 
stars with blue excess and red excess in comparison to a two group population, we find the 
ratio is in fact ~ 1.0 — 1.2, consistent with random errors. Therefore, the assumption that 
there are in fact two distinct stellar populations, K and M giants and supergiants, and OB 
stars, appears to be quantitatively and qualitatively the correct assumption. 

One subject of interest to Galactic Center researchers is whether there is a difference 
between the stellar populations and the extinction (as derived from the stellar population), 
as a function of position on and off the thermal radio emission regions. We explore this 
question, but given the constraints necessary in using the J band data, note that the derived 
ratios should be treated as upper limits. For each of the 70" x 70" images in the El E2 
Filaments and QPS, it was determined whether the position of the frame was coincident 



- 16 - 


with predominately thermal radio emission (the Filaments, Pistol or Sickle), non-thermal 
emission (the Straight Filaments) or did not correspond to any radio features. For the 
H1-H8 region, since there is no non-thermal emission, the regions were divided into the 
thermal emission, the southeast corner and the northwest corner. Since there are only 9 
frames for Sgr A East, we looked at the data in columns of ~ lg W x 3^ s . The results of 
this analysis are presented in Table 4, which gives the median H-K' and J-H values, and 
the average number of stars per frame based on whether they are included on the OB or 
KM extinction curve. To check on the significance of any determined variations, we look at 
random spatial variations by performing the same analysis going from north to south in 
rows of ~2' (roughly the distance covered by two images). 

In the QPS region, the extinction between the various regions shows that there is no 
significant increase in extinction on the thermal filaments as opposed to off the filaments, or 
on the straight filaments. Actually, the differences in extinction are more significant simply 
moving from north to south than any possibly induced by location on or off the radio 
regions. In addition, the value of the extinction is far too small for the stars to be interior 
to molecular clouds in the region as Simpson et al. (1997) estimated that the extinction due 
to the molecular cloud is Ay~70. Therefore, we conclude that the extinction to the stars is 
not effected by the thermal or non-thermal Filaments. In the El & E2 filaments, a similar 
conclusion is reached, although more variations in the extinction are seen. In Table 4, we 
see that the largest extinctions are found in those regions which are off any radio feature. In 
addition, we see a definite trend of increasing extinction going from north to south, with an 
overall increase of 5A V ^5. In H1-H8, once again we see little change with location on and 
off the thermal regions, but there is an increase of Ay^4 going from the southeast corner 
to the northwest corner. In Sgr A East, going from east to west, we see a decrease in the 
extinction, with a total maximum change in the extinction of Ay~5. 
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The average number of OB and KM stars per frame (with each frame being 
approximately a square arcminute), is also presented in Table 4. These results are more 
interesting, but given the small numbers of stars which fit all the imposed constraints, they 
should be treated with extreme caution. The most interesting result is the total ratio of OB 
to KM stars in the different regions. From a compilation of studies of the stellar populations 
in the solar neighborhood, Robin & Creze (1986) found that to a limiting magnitude of 
My=-0.5, which corresponds to our limiting J magnitudes for a minimum of 20 magnitudes 
of visual extinction and the intrinsic colors of the stars under consideration, the ratio 
of OB/KM stars is ~ 1.0 ± 0.05. In the QPS region we find OB/KM~0.66 (estimated 
errors are ±0.05 for all regions), in the El & E2 Filaments OB/KM~0.53, in H1-H8, 
OB/KM~0.51; but in Sgr A East, OB/KM~0.25, although the number of KM stars in Sgr 
A East is consistent with the other regions. In QPS, a region we know to have numerous 
OB stars from spectroscopic observations (Figer et al. 1996) in do we see an increase in the 
ratio of OB/KM stars, which indicates that although limited, this method does confirm 
previously known results. In addition, in all of the regions which contain a significant 
amount of thermal radio emission, we see a significant OB population relative to the KM 
population. Although none of the regions have the same ratio as the solar neighborhood, 
there is a significant number of OB stars that must be taken into account when determining 
the intrinsic color of the overall stellar population. In the QPS region, a few additional 
trends are noteworthy; the ratio of OB/KM stars coincident with the straight filaments is 
~0.4, significantly lower than the other regions. On the other hand, in regions where there 
is no corresponding radio emission, OB/KM^,1.0, which is consistent with the findings of 
Figer et al. (1996), as most of the identified hot massive stars he in the region between the 
Pistol and Sickle, off of the primary radio features. 

Finally, we examine the average extinction in Ay from Table 4 in the various regions 
by using the derived H-K' or J-H values, the corresponding stellar type and the following 
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Av = oki [{H ~ K>) ” {H% ~ m 

Av = oIo4 [(J “^ )_(Ji ~ Ft)] 

The constants are derived using the extinction curve of He et al. (1995), adapted to K 7 
using equation (1). Using the same intrinsic colors, ( H — K[) and (J* — H l ), used before for 
Figures 7-8, we then derived the visual extinction. For example, for H-K'~1.59, the median 
value for QPS, on the thermal filaments, for the OB stars we find that Ay~32. In general 
we find that the visual extinction as determined from OB stars exceeds that determined 
by KM stars by an average of Ay^ 6. When combined with the previous analysis which 
indicates the ratio of OB/KM stars is ~0.5 in the regions observed, this results suggests 
that assuming the stars observed in the GC at NIR wavelengths are K and M giants 
and supergiants underestimates the visual extinction by ~ 2 — 3 magnitudes (see §5). In 
converting the average H-K' values to an extinction in A v , one should therefore assume the 
intrinsic magnitudes are a combination of an approximately half as many OB stars as KM 
stars. Finally, the derived extinction from the median J-H, is consistent with that derived 
from the median H-K' to within Ay ^ 0.7. 

3.3. Color-Magnitude Diagrams 

We present the K' vs. J-K' diagrams for all four regions in Figure 9, and the K' vs. 
H-K' diagrams in Figure 10. The Giant-branch curve for the J-K' diagrams is taken from 
Tiede, Frogel & Terndrup (1995) who found that for Baade’s Window, the Giant-branch 
was best fit by: 

(J — K) = (2.00 ± 0.05) - (0.113 ± 0.005)AT (4) 

which was converted to K' using equation (1). Consistent with previous observations of 
the Galactic Center stars, we see that the observed stars scatter along the Giant-branch 
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curve, although at a variety of extinctions. The Giant-branch curve in Figure 10 was taken 
from observations of 47 Tuc extended to higher luminosities by the observations of stars in 
Baade’s window (Catchpole et al. 1990). 

When comparing Figures 9 and 10 a few things are notable. First, we see that the 
average extinction indicated using the J vs. J-K' data is consistently lower that than 
indicated by the K' vs. H-K' diagram. This supports the assumption that the best estimate 
of the extinction to the stars in the field should be derived from the H and K' data only to 
avoid observational bias. Second, in Figure 9, we note there is a paucity of stars between 
Ay~ 10 — 20, especially noticeable in the QPS region where we have a larger sample. This 
suggests we are seeing a foreground population separated from the stars at the Galactic 
Center by intervening dust and gas, consistent with the proposed clouds within 4 kpc of 
earth (Sanders, Scoville & Solomon 1985). 


4. Spectral Imaging 

The spectral images extracted from the data cubes provide two important pieces of 
information: the location of emission line stars and an estimate of the extinction towards 
the sources in question. Three color images of using emission images derived from the data 
cubes are presented in Figures 11- 18. In the Figures Br 7 is blue, K' is red, H is green, and 
the contours are the 6 cm radio maps. The K' and H band images are taken from the large 
mosaic images and aligned with 2.2 /im continuum images extracted from the cube. Several 
of the emission line images show the location of hot stars, which appear as strong compact 
sources. The nature of these stars and their spectra were presented in Papers I and II. 

The areas which were selected for observations using the wider slit, drift scanning 
method are regions where there is strong radio flux consistent with thermal emission as 
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seen in the 6 cm radio maps (Yusef-Zadeh & Morris 1987a, b; Morris & Yusef-Zadeh 1989). 
Recombination theory predicts that there should be detectable Br 7 2.16 pm emission in 
these regions. In addition, if hot young stars are located within these regions, their strong 
NIR emission lines are detectable as compact sources within the same data scan. From the 
positions of the various features in the IR continuum, radio and Br 7 images, we are able to 
determine if there are exciting stars near or coincident with the radio emission regions, and 
whether or not such stars have extinction values consistent with that of the diffuse gas. 

In order to estimate the extinction to these regions, we must first determine the 
theoretical flux at Br 7 (2.16 pm), f ^ >retlcal , From Rubin (1968), we derive: 

F^or^ai = 8 9 2 6 P a B 0!^ T e 0 - 35 erg cm" 2 s" 1 (5) 

where /? is ratio of Br 7 photons to Lyman continuum photons, a# is the Case B 
recombination coefficient, S„ is the measured radio flux at frequency u, and T e is the 
electron temperature. j3 is calculated using the results for Case B recombination from 
Storey & Hummer (1995) and Hummer & Storey (1987) for the assumed electron density 
and temperature in the Galactic Center (see below). From Rubin (1968) we have that 
o:b=4.10x 10~ 10 Tj 0 - 8 which agrees at T e = 10 4 with the more recent results of Osterbrock 
(1989). The above equation is derived using the same approach as that of Doherty et al. 
(1994) for observations of ultracompact H II regions, and also produces results within 10% 
of those derived using the equation of Smith et al. (1995) for an electron density, n e , equal 
to 10 4 cm -3 and T e = 10 4 K. 

For the Galactic Center, typical values of n e and T e for the known thermal emission 
radio regions vary from region to region. For the Pistol and Sickle, typical n e are ~1000-3000 
cm -3 (Simpson et al.. 1997); for the regions H1-H8, n e range from ~400-3000 cm -3 (Zhao 
et al.. 1993). Since the calculated value of /3 is very weakly dependent upon electron 
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density (varying ^2% for order of magnitude changes in the electron density), we will set 
n e =1000 cm -3 . There is also a range of calculated electron temperatures varying from 5500 
K for the region GO. 095+0. 01 (Erickson et ah. 1991), to 6000 K for the Pistol and Sickle 
(Simpson et ah. 1997), to 4200 K for the region H2 (Zhao et ah. 1993). The calculated 
value of /? is more sensitive to temperature than density, varying by ~25% with a change 
in temperature of 2500 K, we therefore assume T e =6000 K, to be consistent with previous 
GC calculations (see e.g. Morris & Yusef-Zadeh 1989; Simpson et ah 1997). We therefore 
find the best estimate for (3 is 0.0170+0.004 

Finally, for all of the subsequent calculations we use the 6 cm radio maps. Therefore, 
to convert from measured radio flux to predicted Br 7 flux for our Galactic Center data, we 
find from equation (5) that: 

F .theoretical = j 42 x 10~ 21 W cm” 2 (6) 

Using this equation, the extinction is then given by: 

observed ^theoretical —TBr-* (y\ 

r B n ~ r Bn e ' ' 

To determine Fg r s 7 erved , we summed over those regions of Br 7 emission which are 
described in detail in the following sections. The results are presented in Table 5. Rather 
than rely on published flux values, we obtained the copies of the 6 cm VLA radio maps 
available at http://imagelib.ncsa.uiuc.edu/project/document/95.FY, in order to insure the 
closest correspondence possible between the areas summed over in both the radio and our 
Br 7 images. Comparison of our derived values with those published for H5(W), HI and H2 
by Zhao et al. (1993), indicate our derived radio fluxes are correct to within 10%, which 
results in an extinction error of ±0.1 dex. The derived values for rg r7 are presented in 
Table 5. These values of tb t1 give the extinction directly to the radio emission regions, with 
no assumptions made as to galactic location or extinction curves. As such they provide a 
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completely independent estimate of the extinction to the radio sources, which are believed 
to be at the GC (e.g. Yusef-Zadeh & Morris 1987a). Comparisons with the extinction 
derived from the stellar population and discussion of the individual sources are presented 
below. 

In addition to Br 7 , we detected He I (2.058 /jm) emission in almost every region. 
However, recent work by Shields (1993) and DePoy & Shields (1994) on the use of the 
He I/Br 7 ratio as a nebular diagnostic indicates that a single ratio value is found for 
significantly different nebular parameters; hence, unless the electron temperature and 
density are well known, use of this ratio to determine the temperature of the ionizing 
source was not particularly effective. In addition, there are problems properly subtracting 
the atmosphere at the 2.058 fj , m line, especially during long drift scans, thereby further 
reducing the quality of any derived ratio. Finally, Shields (1993) found that the dust 
content and internal velocity distributions significantly effected the line ratio. Therefore, 
since the requisite detailed information is not known for the regions under consideration, 
and we are dealing with complex morphologies, rather than the more simple planetary 
nebulae investigated in DePoy & Shields (1994), the He I/Br 7 is not used in the analysis of 
these regions. 


4.1. Quintuplet, Pistol and Sickle 

Three data scans were taken in the QPS region; the one taken at the southern end 
of the Sickle did not show Br 7 emission significantly above our noise and therefore is not 
presented. In addition to the emission from the Pistol itself, the emission image of the 
Pistol revealed four emission line stars discussed in detail in Paper I. In the combined Bry 
H, K', and radio image, presented in Figure 11, we clearly see that the diffuse emission 
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from the Pistol closely follows the radio emission. Directly above the Pistol, there are 
two compact sources (numbers 8 and 10) w r hich correspond to two stars that are part of a 
group of stars known alternately as AFGL 2004, or the Quintuplet Cluster after the five 
luminous cool stars (numbered 1, 2, 3, 4 and 9 in Figure 11), originally detected in K band 
images (e.g. Nagata et al. 1990). In addition to #8 and #10 directly above the Pistol, 
we see a compact source in the extreme upper left corner corner, corresponding to a star, 
which we call Source B. Directly below the Pistol we see a compact source, we call Source 
A in Paper I, which has recently been observed with HST (Figer et al. 1998), and is now 
known as the ‘Pistol Star’. The Pistol Star is a single star which may be one of the most 
massive stars yet detected (Figer et al. 1998), and corresponds to a slight peak in the radio 
emission (Paper I). In comparison to the other stellar sources in our image it appears to be 
slightly extended. The spectra and discussion of the stars is presented in Paper I and is 
also discussed by Figer et al. (1996) and Morris et al. (1996). We also detected significant 
He I 2.058 /urn emission from the Pistol which follows the Br 7 emission. Notable in the He I 
image, however, is the lack of emission from the Pistol Star. 

The Br 7 emission image of the Sickle (GO. 18-0.04, Figure 12), is interesting for a 
number of reasons. As seen in the other images, the morphology of the Bry emission is 
similar to the radio emission, both of which are extremely clumpy. Although the emission 
comes from a series of compact knots, none of the knots correspond to stars in the 
continuum images. In addition, there are no emission stars detected in this scan although 
there may be additional OB stars nearby (see Figer et al. 1996). The Sickle is the only 
scan in this data set to cover portions of the Straight Filaments (see Figure 1) which are 
approximately constant in declination at scale, so would be parallel to the x-axis in the 
image. No emission from the Straight Filaments is seen, consistent with the nonthermal 
nature of that emission. 



- 24 - 


4.2. El and E2 Filaments 

The radio emission region GO. 10+0.02 is the site of the strongest 6 cm radio emission 
in the Arched Filaments and the radio map also shows several clumps of concentrated flux 
within the region. From the observed FIR emission lines Erickson et al. (1991) estimated 
that numerous O stars are needed to ionize the region. The Br 7 emission image of the 
region combined with the H and K' images and overlaid with the radio contours is shown 
in Figure 13. In the image we see that the diffuse emission is again similar to the radio 
emission. In the upper left region of the image, there is a compact Br 7 emission source, 
labeled ‘A which is coincident with a star . The spectrum of the star is discussed in 
detail in Paper II and appears to be either an extreme Of star or a late type WN star. As 
can be seen in Figure 13, however, no emission line stars are embedded within the radio 
region itself. The spectra of several of the stars lying interior to the bottom of GO. 10+0.02 
were extracted from the data cube of the region: none exhibit emission lines, and all had 
indications of CO absorption at ^2.3 fi m, indicating they are cooler M giants and therefore 
unlikely to be the ionization sources for the region. 

Only one region north of GO. 10+0.02 along the E2 Filament was observed with a 
spectral scan. The location of the scan, labeled “E2 Filament” in Figure 1, is at a peak 
in the 6 cm radio emission. The three color Br 7 emission image is presented in Figure 14. 
Once again the Br 7 emission has the same appearance as the radio emission; however, in 
the E2 Filament Br 7 image, no compact emission sources are seen. In both GO. 10+0.02 and 
the E2 filaments, the He I emission images at 2.06 //m were extracted from their respective 
data cubes but yielded little additional information as the flux was not significantly above 
the background. 
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4.3. Regions H1-H8 

H1-H8, a group of compact H II regions between Sgr A and the Arched Filaments, 
do not appear to be morphologically related to these regions other than by their apparent 
proximity. To date, little IR work has been done on the H1-H8, although they are believed 
to be in the GC (Yusef-Zadeh & Morris 1987b; Zhao et al. 1993). In HI, presented in Figure 
15, we see once again see the similarity of the diffuse emission and the radio emission. 
Although there are stars within the region, none show significant line emission above the 
diffuse background level. There are stars located coincidentally along the arc of diffuse 
emission, but when compared to the stars found in the other images do not stand out. This 
is the only region in this study where the emission mapped in Br 7 and the radio is also 
detected in the K 7 image. In the extreme right of Figure 15, the Br 7 emission corresponds 
to region H2 which is shown more clearly in Figure 16. In the Br 7 image of the H2 region, 
the radio morphology is again reproduced; this is particularly apparent in the arc feature 
below the bright compact region which seems to connect to the H2 region. There are two 
compact emission features at the peak of the H2 emission, one of which is likely a B[e] 
star and the other a Red Giant or Supergiant, both discussed in Paper II. At the peak of 
H2, neither the Br 7 or the radio emission are perfectly spherical, but instead are slightly 
elliptical, with the major axis of the ellipse laying on an axis between the two stars (as can 
clearly be seen in Figure 16). Therefore, it seems likely that both stars are associated with 
the region H2. 

Nearer to the Arched Filaments, observations were made of the region including the 
western section of H5, and all of H6 and H7. In the Br 7 image (Figure 17), we see an 
interesting emission line star. The star, labeled ‘A ’, is located to the east of the crescent 
of emission from H5 and is nearly at the center of curvature for the entire H5 region (see 
Paper II, Figure 7). In Paper II we argue that based on both morphological and empirical 
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grounds, this star is responsible for the ionization of H5. Also seen in Figure 17 is the 
emission from the nearby regions H6 and H7. The emission from H6 is extremely faint while 
the emission from H7 is more easily seen. It is most likely that H7 is a compact H II region, 
although in the continuum images there are two possible stars which may be associated 
with H7. 

Finally, a drift scan was made of the region including H8 (Figure 18), and a star was 
determined to be located within the region. The spectrum of this star indicates it is very 
likely a B[e] star, a more detailed discussion of its stellar type is presented in Paper II. 


4.4. Sgr A East, Regions A-D 

The radio emission to the east of Sgr A East, regions A-D (Figure 19), comes from 
sources which are much more compact and spherically symmetric than most of the previous 
regions. However, the regions appear to be slightly distorted in the direction of the Sgr 
A East. Yusef-Zadeh & Morris (1987a) have suggested the curvature in the flux is due to 
interaction with the expanding SNR in Sgr A East. In the Br 7 emission image (Figure 19), 
we see both the emission from the probable H II regions, and a compact emission line star 
adjacent to region A. This star is fascinating as it contains some of the strongest, broadest 
lines of any star seen in the area and is almost certainly a WN6 star (see Paper II). To the 
left (east) of region D, we see two faint emission sources, which could be additional stars, 
however, results from the other images suggest that if the emission were stellar, it would 
be much stronger as compared to the diffuse emission. Therefore this faint emission may 
be the result of difficulties in background subtraction as the sky conditions varied over the 
1.5 hour scan integration time. Finally, the regions A-D have the strongest He I 2.058 /zir i 
emission of any of the sources. 
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5. Extinction 

As seen in §3.3, the inclusion of J band data introduces a bias towards foreground stars, 
therefore we use only the H and K' data to derive the extinction to the regions which we 
will compare with the Br7 extinction. In addition, since we will be using the average value 
of H-K', we further restrict the sample to stars with 0.8^H-K'<;3.0, to eliminate obvious 
foreground stars and to avoid the problem of incompleteness as we move to fainter values 
of H. We also do not take stars with H below the limiting magnitude for the given regions 
as presented in Table 1. The average color values by frame are presented Figures 20-23 
overlaid on the radio contours. The error on the determined average color is ~0.02 dex. 

The previous study of this kind in this region (Catchpole et al. 1990), had limiting 
magnitudes of H=12.8 and K=12.1 and assumed the stellar population was comprised 
only of K and M stars. With our limiting magnitudes of K'~14-15 and H ~15-16, we 
are observing 4-18 times more stars than the previous study. If, using our data sets, we 
recalculate the average colors imposing the limiting H magnitude of Catchpole et al.(12.8) 
we find that the derived extinction values decrease by Ay^l-4 magnitudes (assuming a 
population of K and M stars only). Thus, if we restrict our data to the same population 
observed by Catchpole et al. we find the calculated average extinction values decreases by 
15%. However, we have argued (§3.2) that the stellar population is a mixture of OB and 
KM stars. Assuming a mixed population we find the visual extinction given by Catchpole 
et al. is underestimated by 4-7 magnitudes, or by ~20%. 

In §3.2, we saw that, except for the region near Sgr A East, the stellar population is 
divided approximately into twice as many K and M giants and supergiants as OB stars. 
Therefore, to derive the extinction in Ak from the numbers given in Figures 20-22, we use 
a weighted average value for these stellar types, ( H x — K[) = 0.11 ± 0.01, with K converted 
K' using equation (1). The derived values of Ak for the regions studied using the Bry 
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maps are presented in Table 6. For the region Sgr A East, since the ratio of OB/KM stars 
~0.25, and we have very few stars on which to make this assessment, we have assumed 
that the stars are best represented by K and M giants and supergiants, and therefore use 
(Hj — K() km — 0.16 to calculate Aj< for regions A-D. 

As can be seen in Table 6 the values derived from the Br 7 studies, except in the regions 
HI and H6, are consistenly lower than those derived by looking at the corresponding stellar 
populations, up to a factor of two in a few cases. One possible explanation for this is that 
the values of the radio flux taken from the VLA map are underestimated due to the losses 
inherent in interferometric observations. To check on this possibility, we have used the single 
dish 43 GHz map of Sofue et al. (1986) to derive the radio flux to the two regions in this 
study which can be isolated in the 43 GHz map: GO. 10+0.02 and the Sickle. We find that 
for GO. 10+0.02, from the 43 GHz map r^ n ~ 2.6 as compared to 2.1 from the 6 cm radio 
map. For the Sickle, we find that t& n derived from the 43 GHz map is the same as that 
derived from the 6 cm map. These comparisons, however, suffer from differences in beam 
size and an inability to accurately match corresponding regions in our NIR observations 
and the 43 GHZ observations. Therefore, as a second check, we calculated the extinction 
we would derive from our Br 7 data if 50% of the VLA emission was lost. The results of this 
comparison are listed in column 3 of Table 6. By assuming that there is emission missing in 
the VLA radio observations, we can reconcile our estimates of the extinction derived from 
the measured Bry emission to that derived from the stellar population in the Pistol, Sickle, 
HI, H5-H7, and D. 

The derived extinctions in regions H2 and H8, however, cannot be reconciled by 
assuming radio losses. This is most likely a result of the Br 7 emission in these two 
regions having a strong stellar component. In these cases, it is not expected that Case B 
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recombination theory would provide an accurate estimate of the Bry emission (Thompson 
1984). The fact that the extinction values in region D can be reconciled assuming only radio 
losses supports the idea that the emission lines seen coincident with a star at that location 
(Paper II) are nebular rather than stellar in origin since a strong stellar contribution would 
suggest a greater discrepancy in extinction estimates. 

The remaining regions, GO. 10+0.02, the peak of the E2 Filament, and A-C, are more 
difficult to understand. One possibility is that these regions are foreground objects. For 
GO. 10+0.02 and the E2 peak, however, this is not a particularly good explanation as 
the Arched Filaments are apparently interacting with the Straight Filaments (Morris & 
Yusef-Zadeh 1989), and the observed radio emission is likley on the surface of a dense 
molecular cloud which we would not anticipate penetrating (Serabyn &: Giisten 1987). 

In addition, the data from GO. 10+002 and the peak of the E2 Filament have the lowest 
signal-to-noise data in the set, and the derived analytical errors may have underestimated 
the quality of the measurements. In the regions A-C, however, no such problem exists, and 
the possibility that these regions are in fact foreground objects cannot be ruled out by this 
analysis. 

The primary conclusion that can be drawn from this comparison is that the ionized gas 
seen in the radio maps cannot be embedded within the coincident molecular clouds. The 
extinction within the +25 km s -1 cloud, co-located with the Sickle, has been estimated to 
be Ay~70 (Simpson et al. 1997), therefore no stars would be observable behind this cloud. 
In addition, if the ionized gas observed in the radio and IR were interior to this cloud, 
the extinction derived from the Br7 data, would be expected to be much larger, but again 
this is not the case. The observed ionized gas is therefore most likely on the front face of 
the molecular clouds as suggested earlier by Serabyn & Giisten (1987). However, with the 
effects from possible underestimates of the VLA radio flux and the difficulties in matching 
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the areas observed with the 43 GHz map with the data presented here, determining whether 
the ionized gas is also in front of part of the stellar populations is problematic at best. 
Finally, given the potential errors in the extinction as estimated from a comparison of radio 
and Br 7 data discussed above, the extinction value we derive from the stellar population 
appears to be the more consistent and reliable of the two methods. 


6. Conclusions 

There appear to be no correlations between the variations in the local stellar population 
and the location of the radio emission regions: there are no increases in stellar density with 
location on or off the radio regions, and no apparent effect in the average stellar extinction 
with location on the thermal emission regions. Studies of the stellar population indicate 
that the ratio of OB to KM stars is about half that of the solar neighborhood in several 
regions. This contradicts the assumption that all of the stars in the GC are evolved K 
and M giants and supergiants, and provides further evidence that there has been recent 
star formation in this region. In addition, previous extinction estimates may be too low 
by 20% due to the limiting H and K magnitudes and the assumed stellar population. 
The derived extinction directly to the radio emission regions is in general, less than the 
stellar population indicating, at a minimum, that the emission seen from these regions in 
on the face of the observed molecular clouds, and some of the regions may be closer than 
the Galactic Center. In addition, the ionization of these regions is stellar in origin as many 
of the regions have either nearby, or embedded emission line stars. 
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Fig. 1. — Position of the infrared images presented in this paper overlaid with the 6 cm radio 
contours (Yusef-Zadeh & Morris 1987a). The shaded areas indicate the regions where large 
broad band images were obtained. The labeled regions indicate radio regions also observed 
in Br7- 

Fig. 2. — Three color J, H, K' image of the QPS region. J is blue, H is green and K' is red. 
The image is a mosaic of 24 70"x70" images, obtained in a 4x6 array. Seeing for the image 
was ^0.8". 

Fig. 3'. — Three color J, H, K' image of the H1-H8 region. J is blue, H is green and K' is 
red. The image is a mosaic of 30 70"x70" images, obtained in a 5x6 array. 

Fig. 4. — Three color J, H, K' image of the Sgr A East region. J is blue, H is green and K' is 
red. The image is a mosaic of 9 70 ,, x70" images, obtained in a 3x3 array. The enhancement 
of stars seen in the- northwest corner, is the concentration of stars along the plane of the 
galaxy. 

Fig. 5. — The log of the number of stars versus magnitude as detected in each of the four 
regions. As can be seen, many more stars were detected in the QPS region due to the 
subarcsec seeing. The smaller total number of stars in Sgr A East region is due the result of 
having taken fewer images of this region. 

Fig. 6. — Median error in magnitude vs. median K magnitude. The largest errors are in 
the QPS frames, due to undersampling (see text). As expected, the error generally increases 
with increasing magnitude. 

Fig. 7. — The Color-Color diagram for QPS and the El and E2 filaments. Errors bars are 
shown in the lower right corner along with the case used (see text for details and Table 2a 
for numbers). The shaded boxes are the estimated errors applied to the reference stellar 
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extinction curves. Stars within these boxes are assumed to be of the same stellar type as 
the nearest spectral curve. 

Fig. 8. — The Color-Color diagram for the regions H1-H8 and Sgr A East. Errors bars are 
shown in the lower right corner along with the case used (see Table 2b). The shaded boxes 
are the assumed errors applied to the reference stellar extinction curve. 

Fig. 9. — The K' vs. J-K' diagram for all the regions. Errors bars are shown in the lower 
right corner along with the case used (see Tables 2a-2b). The Giant-branch curve is from 
Tiede et al. (1995) and the extinction curve that of He et al. (1995) both reference curves 
have been adjusted to K' using equation (1). 

Fig. 10. — The K' vs. H-K' diagram for all the regions. Errors bars are shown in the lower 
right corner along with the case used (see Tables 2a-2b). The Giant-branch curve is from 
Catchpole et al. (1990) and the extinction curve that of He et al. (1995) both adapted to K' 
using equation (1). 

Fig. 11. — Three color image of GO. 15-0.04, the ‘Pistol’. K' is red, H is green, Bry is blue. 
Overlaid with the VLA 6 cm radio contours (Yusef-Zadeh & Morris 1987a). Numbered stars 
8 and 10, and lettered stars A and B correspond to stars with emission line spectra which 
were discussed in Paper I. 

Fig. 12. — Three color image of GO. 18-0.04, the ‘Sickle’, overlaid with the VLA 6 cm radio 
contours (Yusef-Zadeh & Morris 1987a). K' is red, H is green, Br 7 is blue. 

Fig. 13. — Three color image of GO. 10+0.02. K' is red, H is green, Br 7 is blue. The source 
labeled A in the upper left is a compact source we believe to be an emission line star. The 
spectrum of this hot young star and a discussion of its stellar type is presented in Paper II. 

Fig. 14. — Three color image of the radio peak in the E2 Filament. K' is red, H is green, 
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Br7 is blue. 

Fig. 15. — Three color partial image of the compact regions HI, left, and H2, right. K' is 
red, H is green, Br 7 is blue. See Figure 1 for location relative to the Arched Filaments and 
the Sgr A West halo. No emission line stars are apparent in this image, but the center of 
curvature of HI was not covered by this scan, so the prescence of a hot star in this region is 
not precluded. 

Fig. 16. — Three color image of the compact region H2. K' is red, H is green, Br 7 is blue. 
The star labeled A has an emission spectrum indicative of a B[e] star. The star labeled B 
has the deep CO absorption features suggesting an M supergiant. The stars are discussed 
in more detail in Paper II. 

Fig. 17. — Three color image of the compact regions H5-H7. K' is red, H is green, Br 7 is 
blue. The star located within the emission of H5, labeled A, shows strong compact Br 7 
emission and could be a B[e] or possibly an LBV star. See Paper II for details. 

Fig. 18. — Three color image of the compact region H8. K' is red, H is green, Bry is blue. 
The star seen in conjunction with the region is the only one to shown Bry in emissin and is 
most likely a B[e] star. See Paper II for details. 

Fig. 19. — Three color image of the compact regions A-D located to the east of Sgr A East. 
K' is red, H is green, Br 7 is blue. The emission star marked with the asterisk is a WN 6 star 
as discussed in Paper II. 

Fig. 20. — Average H-K' for the QPS region. Each number is centered on the location of 
the individual image used to derive the color. Each image is ~l'xl'. 

Fig. 21. — Average H-K' for the El & E2 Filaments. Each number is centered on the 
location of the individual image used to derive the color. Each image is ~l , xl'. 
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Fig. 22. — Average H-K' for the H1-H8 region. Each number is centered on the location of 
the individual image used to derive the color. Each image is ~l'xl'. 

Fig. 23. — Average H-K' for the Sgr A East region. Each number is centered on the location 
of the individual image used to derive the color. Each image is ~l'xl'. 








Table 2a. Median Errors for Continuum Images 


QPS El & E2 Filaments 

Number Number 


Case 

Band 

of stars 

Mag 

a 

DAOPHOT 

of stars 

Mag 

a 

DAOPHOT 

2A 

K' 

1659 

13.46 

0.076 

0.065 

1636 

12.00 

0.050 

0.034 


H 


14.95 

0.059 

0.068 


13.81 

0.051 

0.040 

2B 

K' 

515 

12.18 

0.072 

0.040 

303 

10.62 

0.038 

0.018 


J 


16.62 

0.052 

0.061 


15.80 

0.040 

0.040 

2C 

K' 

492 

12.18 

0.071 

0.040 

299 

10.61 

0.040 

0.018 


H 


13.56 

0.047 

0.043 


12.27 

0.037 

0.016 


J 


16.65 

0.052 

0.062 


15.80 

0.038 

0.038 

4A 

K' 

78 

13.02 

0.127 

0.057 

103 

11.81 

0.078 

0.032 


H 


14.65 

0.097 

0.062 


13.43 

0.070 

0.034 

4B 

K' 

24 

11.84 

0.137 

0.040 

18 

10.56 

0.063 

0.017 


J 


16.52 

0.077 

0.060 


15.50 

0.065 

0.030 

4C 

K' 

23 

11.73 

0.130 


17 



0.017 


H 


13.11 

0.082 



11.69 


0.013 


J 


16.47 

0.077 



15.53 


0.029 


Table 2b. Median Errors for Continuum Images 


Case 

Band 

Number 
of stars 

H1-H8 
Mag a 

DAOPHOT 

Number 
of stars 

Sgr A East 
Mag a 

DAOPHOT 

2A 

K' 

1289 

12.87 

0.046 

0.037 

215 

12.18 

0.049 

0.023 


H 


14.65 

0.056 

0.046 


13.90 

0.045 

0.040 

2B 

K' 

316 

11.94 

0.036 

0.026 

62 

11.30 

0.043 

0.016 


J 


16.80 

0.054 

0.043 


16.00 

0.039 

0.030 

2C 

K' 

309 

11.93 

0.035 

0.025 

62 

11.30 

0.043 

0.016 


H 


13.36 

0.036 

0.025 


12.58 

0.033 

0.017 


J 


16.79 

0.052 

0.042 


16.00 

0.039 

0.030 

4A 

K' 

67 

12.86 

0.079 

0.039 

9 

11.45 

0.056 

0.016 


H 


14.57 

0.087 

0.043 


12.91 

0.064 

0.030 

4B 

K' 

13 

11.70 

0.071 

0.030 

1 

. . . 




J 


16.59 

0.075 

0.043 





4C 

K' 

13 

11.71 

0.072 

0.030 

1 





H 


13.36 

0.057 

0.027 






J 


16.58 

0.075 

0.040 








Table 3. Intrinsic Infrared Magnitudes 


Type 

My a 

M k 

M K ' b 

M h 

Mj 


Main Sequence (V) 

05 c 

-5.7 

-4.81 

-4.82 

-4.89 

-4.97 

B0 C 

-4.1 

-3.34 

-3.35 

-3.41 

-3.52 

B5 C 

-1.1 

-0.68 

-0.70 

-0.77 

-0.77 

A0 d 

0.7 

0.70 

0.70 

0.70 

0.70 

F0 d 

2.6 

1.90 

1.91 

1.93 

2.06 

G0 d 

4.4 

2.99 

3.0 

3.04 

3.35 

K0 d 

5.9 

3.94 

3.95 

4.02 

4.47 

M0 d 

9.0 

5.35 

5.37 

5.52 

6.21 

Giants (III) 

05 c 

-5.4 

-4.52 

-4.53 

-4.60 

-4.72 

B5 C 

-5.0 

-4.31 

-4.33 

-4.40 

-4.49 

B0 C 

-2.4 

-1.99 

-2.00 

-2.05 

-2.09 

G0 d 

1.1 

-0.65 

-0.64 

-0.59 

-0.20 

G5 d 

0.7 

-1.40 

-1.39 

-1.32 

-0.84 

K0 d 

0.5 

-1.81 

-1.79 

-1.72 

-1.18 

K5 d 

-0.2 

-3.80 

-3.77 

-3.64 

-2.85 

M0 d 

-0.4 

-4.25 

-4.22 

-4.06 

-3.24 

M5 d 

-0.8 

-6.76 

-6.71 

-6.48 

-5.53 

Supergiants (la) 

05 c 

-8.4 

-7.50 

-7.51 

-7.53 

-7.69 

B0 C 

-6.7 

-6.09 

-6.09 

-6.11 

-6.22 

B5 C 

-6.9 

-6.69 

-6.69 

-6.69 

-6.76 

K5(Ib) e 

-4.6 

-8.10 

8.08 

-7.96 

-7.35 

M0 e 

-7.5 

-11.3 

-11.3 

-11.13 

-10.58 

M2 e 

-7 

-11.1 

-11.1 

-10.89 

-10.34 


a Allen, Astrophysical Quantities 
b Derived from equation (1) 
c Infrared colors from Wegner & Walter (1994) 
d Infrared colors from Bessel & Brett (1988) 
e Infrared colors from Elias et al. (1985) 



Table 4. H-K', J-H, and Average Number of stars by Location and Stellar 

Type 


H-K' J-H Ave. #of stars/Frame 

OB KM OB KM OB KM 


QPS 


On 

1.59±0.01 

1.48±0.01 | 

3.26±0.01 

3.36±0.01 

20 

32 

Off 

1.60±0.02 

1.49±0.01 ! 

3.27±0.05 

3.50±0.01 

21 

20 

Straight 

1.55±0.02 

1.45±0.01 

3.19±0.04 

3.374=0.03 

11 

30 

Row 1 

1.53=b0.01 

1.454:0.01 

3.28±0.03 

3.38±0.02 

14 

22 

Row 2 

1.60±0.01 

1.524:0.01 

3.34±0.03 

3.55±0.01 

23 

29 

Row 3 

1.58±0.01 

1.454:0.01 

3.31±0.02 

3.38±0.02 

26 

30 


El and E2 Filaments 


On 

Off 

Straight 

1.89±0.04 

1.96±0.09 

1.67±0.01 

1.73±0.04 

1.80±0.03 

1.56±0.01 

3.91±0.10 

4.03d=0.19 

3.46±0.01 

3.88±0.07 

4.03±0.05 

3.59±0.01 

5 

6 
8 

10 

12 

15 

Row 1 

1.66±0.01 

1.54±0.02 

3.44±0.01 

3.58±0.01 

8 

16 

Row 2 

1.85±0.04 

1.74±0.04 

3.90±0.14 

3.92±0.08 

9 

13 

Row 3 

1.80±0.02 

1.78±0.03 

3.75±0.01 

3.97±0.04 

4 

12 

Row 4 

1.98±0.10 

1.81±0.10 

4.12±0.23 

4.09±0.22 

5 

7 


H1-H8 


On 

NW 

SE 

1.74±0.02 

1.83±0.02 

1.63±0.05 

1.58±0.01 
1.60±0.02 
1.56±0.01 i 

3.59±0.02 

3.75±0.05 

3.44±0.03 

3.65±0.03 

3.58±0.02 

3.60±0.01 

7 

10 

8 

15 

7 

20 

Row 1 

1.86±0.02 

1.63±0.01 

3.78±0.02 

3.69±0.02 

6 

11 

Row 2 

1.77±0.03 

1.64±0.02 

3.64±0.07 

3r.77±0.07 

8 

14 

Row 3 

1.65±0.04 

1.53±0.02 

3.44±0.05 

3.54±0.01 

9 

22 


Sgr A East 


Column 1 

1.51±0.03 

1.69±0.01 

3.33±0.15 

3.93±0.04 

1 

14 

Column 2 

1.59±0.03 

1.45±0.01 

3.39±0.06 

3.39±0.02 

7 

12 

Column 3 

1.85±0.31 

1.48±0.02 

3.97±0.68 

3.51±0.05 

1 

10 



Table 5. Observed vs. Theoretical Br7 
emission 


Region 

Figure 

Measured Br 7 
(10^ 20 W cm -2 ) 

6 cm 

(mJy) 

TBr'y 




QPS 

Pistol 

mm 

5.73+0.08 

520 

2.6 

Sickle 


12.6±1.2 

1000 

2.4 


! El & E2 Filaments 

GO. 10+0. 02 

■■ 

19.0+0.7 

990 

2.0 

E2 Peak 

m 

19.5±3.0 

260 

0.6 


H1-H8 


HI 

15 

2.2±0.1 

470 

3.4 

H2 

16 

11.7±4 

500 

1.8 

H5(W) 

17 

6.2±0.4 

510 

2.5 

H5(E) 

17 

0.63+0.08 

45 

2.3 

H6 

17 

0.076+0.003 

16 

3.4 

H7 

17 

0.39+0.10 

45 

2.8 

H8 

18 

0.45+0.01 

16 

1.6 


Sgr A East 


A 

19 

17.1±0.4 

420 

1.2 

B 

19 

3.2±0.2 

100 

1.5 

Cl 

19 

3. 1+0.2 

87 

1.4 

C2 

19 

1. 1+0.1 

24 

1.2 

D 

19 

0.63+0.06 

42 

2.3 











Table 6. Comparison of Ak 



1 Photometry 

Derived from 
Br 7 Extinction 

Derived from Br 7 
(assuming radio losses) 

Pistol 

3.0 

2.7 

3.0 

Sickle 

3.0 

2.5 

3.2 

GO. 10+0.02 

4.2 

2.1 

2.8 

E2 Peak 

3.3 

0.6 

1.4 

HI 

3.6 

3.6 

4.3 

H2 

3.7 

1.9 

2.6 

H5 

3.6 

2.5 

3.2 

H6 

3.6 

3.6 

4.3 

H7 

3.6 

2.9 

3.6 

H8 

3.6 

1.7 

2.4 

A 

2.8 

1.3 

2.1 

B 

2.8 

1.6 

2.3 

Cl 

2.8 

1.5 

2.2 

C2 

2.8 

1.3 

1.9 

D 

2.8 

2.4 

3.1 
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